Isospin equilibration in multinucleon transfer reaction at near-barrier energies 
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The isospin equilibration process in multinucleon transfer reactions is investigated by using the 
improved quantum molecular dynamics model. The collisions of !?*Xe4-?95 Pb at near-barrier en- 
ergy are studied with different symmetry energy coefficients. We find that the neutron transfer is 
enhanced at the early stages of the collisions under the strong symmetry potential. The neutrons 
transfer takes place at earlier stage than protons. The trajectory of neutron flow is along the low- 
density path from the target to projectile. The average N/Z values of the primary products in 
58N1.-705pb reaction are also investigated and compared with available experimental data. It shows 
that the average N/Z values of the projectile-like products increase rapidly with increasing mass 
transfer. In the quasi-elastic collisions, the isospin equilibration process is incomplete due to a short 
contact time between the reaction partners. The complete isospin equilibration events take place in 


symmetric quasi-fission reactions. 


I. INTRODUCTION 


'The isospin transport effects in heavy-ion collisions have 
been studied for many years[1-9]. It includes two differ- 
ent facets, named isospin diffusion and drift. The isospin 
diffusion is related to the isospin asymmetry of a system 
in which the projectile and target have different N/Z val- 
ues [10]. The isospin drift is related to the density gra- 
dient which is expected to exist in the low-density neck 
region, even between two identical nuclei [11]. In the 
nuclear reactions, the isospin transport is initiated and 
continues until the system disintegrates or the chemical 
potentials for neutrons and protons in both nuclei become 
equal. If the interaction time between the projectile and 
target is long enough, the system will reach the situation 
of isospin equilibration. 

In heavy-ion collisions at intermediate energies, the 
isospin transport effects have been investigated in mid- 
peripheral collisions by analyzing the N/Z values of the 
reconstructed primary quasi-projectile, quasi-target, and 
midvelocity sources. Such as in the 4°Ca+*8Ca reaction 
[3], one sees that the isospin diffusion effects lead to an 
increase of the N/Z of the quasi-projectile source. In 
the collisions of two identical N/Z values nuclei (such as 
55Ni4-55Ni [11]), it allows us to observe the isospin drift 
phenomenon without contamination from the midveloc- 
ity source. 

At low incident energies, the multinucleon transfer 
(MNT) reaction is one of the most important mecha- 
nism, which has attracted widespread interests in recent 
years both experimentally [12-22] and theoretically [23- 
35]. Unlike at the intermediate incident energies, the 
collisions between the reaction partners are not violent. 
'The major products in MNT reactions only contain the 
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quasi-projectlie and the quasi-target fragments with the 
number of nucleon transfers up to 30 or more. No midve- 
locity fragments are formed. In general, the MNT reac- 
tions include the quasi-elastic, the deep-inelastic, and the 
quasi-fission reactions. These reaction mechanisms can 
be distinguished by analyzing the total-kinetic-energy- 
mass (TKE-Mass) distributions combined with the total 
kinetic energy loss (TKEL) distributions of the products 
[36-38]. The contact times between the reaction partners 
for the quasi-elastic, the deep-inelastic, and the quasi- 
fission reactions are quite different. The MNT reactions 
provide an opportunity to study the isospin equilibra- 
tion process under the different time scales. In the quasi- 
fission reactions, the typical time scales can reach 107?! 
s order of magnitude or longer [39, 40]. In this case, the 
nucleon transfer process may lead to a uniform distribu- 
tion of the N/Z ratio. 

Many experiments about the isospin equilibration ef- 
fects at low-energies have been performed by analyzing 
the N/Z values of reconstructed primary products during 
the last decades [18, 41-43]. For example, the reactions 
of §4Ni+1°Te at Ejay = 275 MeV and *8Ni+?°8Pb at 
Fiay = 345 MeV have been investigated by Królas et 
al. at Laboratori Nazionali di Legnaro (INFN) [18]. It 
is found that the isospin equilibration process is closely 
related to the number of nucleon transfer between the 
reaction partners. Actually, it is difficult to divide the 
contributions of the isospin diffusion and drift at low- 
energies separately. The isospin equilibration is one im- 
portant mechanism in the MNT reactions, which is one 
promising approach to produce new neutron-rich nuclei. 
For example, the new isotopes ?^Ti, 96V, 55,5990, ©'Mn, 
and 995:915Fe were produced by Guerreau et al. through 
a 340 MeV *? Ar beam accelerated by the Orsay ALICE 
accelerator facility bombarding on ??*U target [44]. 

The neck is an important characteristic in MNT re- 
actions, even in the early stage of the fusion reaction. 
It is formed in the dinuclear system and is character- 
ized by sub-saturation densities. The isospin transport 


and the energy dissipation processes of the system are 
closely related to the nucleon transfer between the re- 
action partners through the neck. The dynamical cal- 
culation is required to understand some details of these 
processes. A microscopic dynamics model, the improved 
quantum molecular dynamics (ImQMD) model [45, 46], 
describe the nuclear reactions based on effective nucleon- 
nucleon interaction. The ImQMD model is self-consistent 
for describing the behavior of neck evolution and nucleon 
transport during the collisions. In this work, we apply the 
ImQMD model to investigate the isospin equilibration 
process in the reactions of 1?4Xe+208Pb and °8Ni+?°8Pb 
at near-barrier energies. 

'The structure of this paper is as follows. In Sec. II, 
we briefly introduce the ImQMD model. The results and 
discussion are presented in Sec. III. Finally, the conclu- 
sion is given in Sec. IV. 


II. THE MODEL 


'The ImQMD model is an improved version of the quan- 
tum molecular dynamics (QMD) model [47]. The stan- 
dard Skyrme interaction with the omission of the spin- 
orbit term is adopted for describing the bulk and sur- 
face properties of nuclei [45]. The stochastic two-body 
collision process is added to the time evolution by the 
Hamilton equation of motion. To describe the fermionic 
nature of the N-body system, the Fermi constraint pro- 
posed by Papa et al. in the constrained molecular dy- 
namics (CoMD) model [48, 49] is introduced. It improves 
greatly the stability of an individual nucleus. The final 
state of the two-body collision process is checked so that 
it obeys the Pauli principle. The detailed descriptions of 
the ImQMD model can be found in Refs. [45, 46]. The 
nuclear interaction potential can be read as 
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Here p = pn + pp is the nucleons density. 6 = (Pn — 
Pp)/(Pn + pp) is the isospin asymmetry. The first three 
terms in above expression are obtained from the potential 
energy functional of Skyrme interaction directly. The 
fourth term is the symmetry potential energy including 
the bulk and the surface symmetry potential energies. 
The surface symmetry potential energy term is related to 
the density gradient which is important for describing the 
neutron skin of nuclei. The last term is a small correction 
term. The density distribution function p of a system can 
be read 
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]. (2) 


c» is the wave-packet width of the nucleon in coordinate 
space. The IQ2 parameter sets (see Table 1) adopted in 


this work are the same as in Refs. [24, 32]. The incom- 
pressibility coefficient, Kə, is 195 MeV. This parameter 
sets have been successfully applied on the heavy-ion col- 
lisions in fusion reactions [46, 50], multinucleon transfer 
reactions [24, 32], and ternary breakup reactions|[51]. 


TABLE I. The model parameters (IQ2) adopted in this work. 


a B — Beur & n Cs ks — po 
(MeV) (MeV) (MeV fm?) (MeV) (MeV) (fm?) (fm?) 


—356 303 7/6 7.0 12.5 2/3 32.0 0.08 0.165 


In this work, we set z-axis as the beam direction and 
z-axis as the impact parameter direction. We set the 
wave-packet width o, = 1.3 fm to calculate the isospin 
equilibration process of !?4Xe--209pp and °8Ni+?°SPb. 
The initial distance of the center of mass between the 
projectile and target is 30 fm. 


III. RESULTS AND DISCUSSION 


We first test the ImQMD model for the description 
of the MNT reactions. The collisions of ?5Ni--?995Pb at 
Ep = 328.4 MeV are simulated by using the ImQMD 
model. 39000 simulation events are calculated with the 
impact parameters from 0 to ba; fm. bmax = Rp + Rr, 
where Rp and Rr denote the radius of the projectile and 
the target, respectively. For each event, we simulate the 
collision process until t = 2000 fm/c with a step size of 
At = 1fm/c. Figure 1 shows the isotope production cross 
sections from Mn to Ni. The fusion and elastic scattering 
events are not taken into account. The experimental data 
are taken from Ref. [52]. The thick folding lines and thick 
solid lines denote the calculations from the combination 
of ImQMD--GEMINI and GRAZING model [53-55] with 
inclusion of the evaporation, respectively. The nuclear 
level densities in the GEMINI code [56] are taken as a 
Fermi-gas form with the default parameters. From Fig. 
1 we find that the measured isotope distributions for the 
elements Co and Ni can be reasonably well reproduced 
by the ImQMD+GEMINI calculations. The discrepancy 
between the data and the calculations increases with the 
increasing number of proton transfer. The GRAZING 
calculations also show a similar behavior. 

For simplicity, we consider head-on collisions of 
124X%e428Pb at Eem. = 450 MeV to investigate the 
isospin equilibration process. For this reaction, the in- 
cident energy is slightly higher than the Coulomb bar- 
rier. It is almost impossible to form a compound nu- 
cleus, due to large ZpZr value. The typical reaction 
process is that the colliding nuclei exchange a num- 
ber of nucleons and then re-separate. Figure 2 shows 
the single-particle potentials of neutrons and protons in 
124X e--205Pb along the beam direction at t = 200 and 
300 fm/c. The single-particle potential is defined as 
Vlr) = f p(r')V(r — r')dr', where q = n,p, p(r) is 
the density distribution of the system, and V(r — r’) is 
the effective nucleon-nucleon interaction. The Coulomb 
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FIG. 1. (Color online.) Isotope production cross sections from Mn to Ni in the reaction of 58Ni.L?08pb at Flap = 328.4 MeV. 
The thick folding lines and thick solid lines denote the calculations from the combination of mQMD+GEMINI and GRAZING 
model with inclusion of the evaporation, respectively. The thin folding lines and thin dash-dotted lines denote the primary 
fragment distributions from the ImQMD and GRAZING models, respectively. The experimental data (solid circles) are taken 


from Ref. [52]. 


interaction is very important for the isospin transport 
processes in the early stage of the MNT reactions. For 
example, at t — 200 fm/c, the neutron transfer is allowed 
due to lower barrier between the two reaction partners 
(see Fig. 2(a)). However, for protons, the proton transfer 
between two reaction partners is forbidden at this time 
because the barrier is very high (see Fig. 2(c)). Hence, 
the neutrons transfer process takes place at earlier stage 
than the protons in the MNT reactions. At t = 300 fm/c, 
the barriers between the reaction partners are reduced 
both for neutrons and protons. Therefore, both the neu- 
trons and protons are allowed to transfer between the 
target and projectile (see Fig. 2(b) and 2(d)). 

Figure 3 shows the time evolution of the density pro- 
files for !?^Xe--?085ppb at Eom. = 450 MeV. As seen from 
the figure, the central densities of nuclei in the reactions 
are reasonable. In addition, we note that the diffuse- 


ness of density can be seen clearly on the surface of the 
system. In Fig. 4, we show the contour plots of isospin 
asymmetry (5 = (p, — pp)/ (pn ^- pp)) for the !?*Xe--?08 pp 
during the evolution of the system. The average values 
of isospin asymmetry for !?^Xe and ?95Pb nuclei are 0.13 
and 0.21, respectively. From Fig. 4, one sees that the dis- 
tributions of the isospin asymmetry in the cores of two 
reaction partners are uniform. The corresponding values 
of 6 in the quasi-projectile and quasi-target are 0.1 and 
0.17, respectively. However, the values of isospin asym- 
metry on the surface of the nuclei are significantly larger 
than those in cores. The core of ?°°Pb is covered by neu- 
tron skin with maximal 6 = 0.54 at t = 50 fm/c. The 
neutron enrichment can also be seen clearly in the low- 
density neck region. This phenomenon is known as the 
isospin drift which is due to the density gradient. The 
neutrons are preferably driven to the low density area 
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FIG. 2. (Color online.) The single-particle potentials of neutrons (up panels) and protons (bottom panels) in !?*Xe4-?95pp at 


t = 200 and t = 300 fm/c. 


[57]. In Fig. 4, one also can see that the neutron flow in 
the reactions is formed along the low-density path. 

In order to investigate the isospin transport during the 
evolution of the dinuclear system, we introduce the sep- 
aration plane to divide the quasi-projectile and quasi- 
target nuclei. The separation plane can be defined as the 
plane at the position where iso-contours of the projec- 
tile and target densities cross each other. This method 
was adopted by the TDHF calculations in Refs. [58, 59]. 
Figure 5 shows the transfer coefficients of neutron and 
proton in the head-on collisions of !?4^Xe--?085pb with 
the symmetry energy coefficients of C, — 28 and 32 
MeV. The transfer coefficient can be read as Un p = 
dN, ,/dt, where Np p denotes the net neutron or pro- 
ton flux through the separation plane. The direction of 
nucleon transfer is from the target to projectile. One can 
see that the neutron transfer is enhanced in the case of 
C, = 32 MeV at the early stages of the collisions due to 
the strong symmetry potential. The corresponding peak 
value of the neutron transfer coefficient is about 1 x 107? 
per fm/c which is close to the TDHF calculations in Ref. 
[60]. The proton transfer coefficients are not sensitive 
to the symmetry energy in !?^Xe--?95Pb reactions. Be- 
cause the symmetry energy mainly influence the neutron 
transport in the case that projectile and target have large 


difference of N/Z values. In addition, we note that the 
time of neutron transfer starts at t = 150 fm/c which 
is earlier than that of proton transfer. This is due to 
that the proton transfer is prevented by a high barrier 
between the reaction partners at the early stages of the 
collisions (see Fig. 2(a) and 2(c)). 

In MNT reactions, the contact time between the reac- 
tion partners is related to collision mode. Different col- 
lision mechanisms can be distinguished by analyzing the 
TKE-Mass distribution combined with the TKEL distri- 
bution of the products. In Fig. 6, we show the T'KE-Mass 
distributions and corresponding TKEL distributions of 
primary binary fragments in ?9Ni--295Pb at Ej = 345 
MeV for the contact time ranges of 0 < teon < 400, 
400 < teon € 1000, and 0 < teon < 1000 fm/c. The 
contour plots show the counts of the fragments in log- 
arithmic scale. The region of the impact parameters in 
calculations is from 0 to 54,4; fm. The quasi-elastic and 
the deep-inelastic collisions occur when the contact time 
is less than 400 fm/c. As shown in Fig. 6 (b), the quasi- 
elastic collision events can be extracted through fitting 
the peak of TKEL distribution by a Gaussian curve. In 
general, the contact time between the reaction partners 
for the quasi-elastic collision events is shorter than that 
for the deep-inelastic collision events. There are only 
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FIG. 3. (Color online.) Time evolution of the density profiles for !?*Xe--?95Pp at Eom. = 450 MeV. 


a few nucleon transfers between the projectile and tar- 
get. One can see that the TKEL of the quasi-elastic 
collision events is less than 30 MeV. The differences be- 
tween the ImQMD calculations and Gaussian curve are 
mainly the deep-inelastic collision events. Most of quasi- 
fission events occur at 400 < teon < 1000 fm/c. The 
masses distribute in a rather broad range. There are a 
large number of nucleon transfers, in which may leads to 
the isospin equilibration between the quasi-projectile and 
quasi-target fragments. 

The average N/Z values of the primary products in 
58Ni4-208pPb at Elab = 345 MeV are shown in Fig. 7. 
The solid line in Fig. 7(a) denotes the calculations of the 
ImQMD model with 0 < teon < 1000 fm/c in the region 
of the impact parameters from 0 to bmax fm. The elas- 
tic scattering and fusion-fission events are not taken into 
account. The N/Z values of the projectile, the target, 
and the compound nucleus are 1.071, 1.537, and 1.418, 
respectively. From Fig. 7(a), one sees that the exper- 
imental data can be reasonably well reproduced by the 
ImQMD calculations. The average N/Z values of the 
projectile-like products increase rapidly with increasing 
nucleon transfer due to the isospin transport from tar- 
get to projectile. On the curve, a steep valley and a 
peak can be clearly seen in the vicinity of the projec- 
tile and target. This feature was also found in the re- 
actions of 94Ni4-? Te and 94Ni4-295Pp [18, 43]. In fact, 


the isospin equilibration process depends strongly on the 
contact time of the reaction partners. Figure 7(b) shows 
the average N/Z values of the primary products with 
contact time regions of 0 < teon < 200, 200 < toon < 400, 
and 400 < tecon < 1000 fm/c. In case of 0 < teon < 200 
fm/c, most of the primary products are produced in the 
quasi-elastic collisions. The strong neutron flow from the 
target to projectile results in sharp increase of N/Z val- 
ues of the projectile-like fragments. However, the system 
can not reach the situation of complete isospin equili- 
bration due to a short contact time. With increasing 
contact time, the enhanced effect of the isospin equili- 
bration can be seen clearly. The products of complete 
isospin equilibration are produced in symmetric quasi- 
fission reactions. The isospin equilibration process in 
the MNT reactions can result in the production of very 
neutron-rich projectile-like fragments. A strong absorp- 
tion of neutrons by the projectile (?5Ni) was observed by 
the experiment [18]. For example, after neutron evapo- 
ration, 9^Ni with production cross section about 15 ub 
has been detected in the y — y coincidence analysis. 


IV. CONCLUSIONS 


In summary, the production cross sections of isotopes in 
the reaction of 95Ni4-298Pp at Ej = 328.4 MeV are cal- 
culated by the ImQMD model. The results show that 
the ImQMD model is suitable to describe the MNT reac- 
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FIG. 4. (Color online.) Time evolution of the isospin asymmetry 6 = (pn — pp)/(Pn + pp) distributions for 1?*Xe+?°8Pb at 


tions at near-barrier energies. The isospin equilibration 
process in !?4Xe+?°8Pb at head-on collisions is studied. 
'The values of isospin asymmetry in the cores of two re- 
action partners are uniform during the collisions. The 
neutrons are preferably driven to the low-density area. 
On the surface of the nuclei, the value of isospin asym- 
metry is larger than the core of the nuclei. The direction 
of neutron flow is from the target to projectile along the 
low-density path. And the neutrons transfer process is 
very sensitive to the symmetry energy. The neutron flow 
is enhanced under the large symmetry energy coefficient. 
The isospin equilibration process depends strongly on 
the contact time of the reaction partners. By analysing 
average N/Z values of the primary binary products in 
58Nij--205Pb at Elab = 345 MeV, we find that the aver- 
age N/Z values of the projectile-like products increase 
rapidly with increasing nucleon transfer. In the quasi- 
elastic collisions, the system can not reach the situation 
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